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Abstract 

The distribution of the azimuthal angle for the charged hadrons has 
been studied in the hadronic centre-of-mass system for neutral current 
deep inelastic positron-proton scattering with the ZEUS detector at HERA. 
Measurements of the dependence of the moments of this distribution on 
the transverse momenta of the charged hadrons are presented. Asymme- 
tries that can be unambiguously attributed to perturbative QCD processes 
have been observed for the first time. 
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1 Introduction 



Semi-inclusive processes in deep inelastic scattering (DIS) are of importance 
because they can be used to test the perturbative Quantum Chromodynamic 
(QCD) description of hadron production via parton fragmentation. An observ- 
able of particular interest is the distribution of the azimuthal angle, </>, (mea- 
sured in the hadronic centre-of-mass frame, HCM) between the lepton scat- 
tering plane, defined by the incoming and outgoing lepton momenta, and the 
hadron production plane, defined by the exchanged virtual boson and an out- 
going hadron (Fig. |l|). 

Asymmetries in the 4> distribution, i.e. terms proportional to cos0 and cos2</> 
arise whenever a non-zero transverse momentum, in the HCM frame, is present 
in the scattering process. Consequently both non-perturbative and perturbative 
QCD effects | | give rise to these asymmetries. The azimuthal dependence of 
parton production has the form 



da/d<j) = A + B cos <p + C cos 2</>. (1) 

This form results from the polarisation of the exchanged virtual boson. The 
coefficients B and C depend on the hclicitics of the final-state parton(s), par- 
tonic transverse momenta and on colour coherence J5| . The cos 20 term is ex- 
pected from interference of amplitudes arising from the +1 and —1 helicity 
components of the transversely-polarised part of the exchanged boson, whereas 
transverse/longitudinal interference gives rise to the cos0 term. 

The asymmetry that results from the intrinsic momentum of a quark in the 
proton is referred to as the non-perturbative asymmetry. Since the intrinsic mo- 
mentum is small, this asymmetry should fall rapidly with increasing transverse 
momentum of the measured hadron and with increasing Q 2 [Q] , where Q 2 = —q 2 
is the negative square of the four-momentum of the virtual exchanged boson. 
In particular, the term C is small at high Q 2 for the non-perturbative effect. 

In contrast, the asymmetry associated with leading-order terms in perturbative 
QCD calculations, the perturbative asymmetry, is weakly dependent on Q 2 
and persists at high transverse momenta. The perturbative QCD contribution 
to terms B and C is large at leading order (LO) in a s . At this order, two 
processes contribute to DIS: QCD-Compton scattering (QCDC, j*q — > qg) and 
boson-gluon fusion (BGF, j*g — > qq). For the QCDC process the scattered 
quark preferentially appears at 4> near 180° whilst for the BGF process the </> 
dependence has a symmetry about 180°, as a consequence of the quark/anti- 
quark final state. The BGF process is the dominant contribution to the cos2</> 
term, while the QCDC contribution dominates the cos <fi term. The cos 2</> term 
can be unambiguously attributed to perturbative QCD processes in the high Q 2 
region under investigation. 

Terms in sin </> and sin 2cf> are also theoretically expected only for deep inelastic 
scattering with polarised leptons or charged current reactions. Even in these 
cases, their magnitude is however estimated to be substantially smaller than 
terms in cos </> and cos 2(f> . 

The form of Eq. |l| is expected to be maintained for single particle produc- 
tion j| , since high- momentum hadrons are produced close to the direction of 
the parton. Measurable perturbative and non-perturbative predictions for the 
asymmetries in hadron production require the use of an appropriate fragmen- 
tation function or hadronisation mechanism. Note that in order to observe the 
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cos 4> asymmetry, a selection procedure is required which consistently associates 
leading hadrons produced from either quarks or gluons |^,|7). This is described 
in Section ^|. The transverse momenta arising from the fragmentation itself does 
not contribute to the asymmetry but only smears the observed distribution Q . 

The event kinematics of DIS are determined by Q 2 , and one of the two Bjorken 
scaling variables x — Q 2 /2P-q or y = Q 2 /xs, where P is the four-momentum of 
the incoming proton and t/s is the positron-proton centre-of-mass energy. The 
kinematic region studied is 0.2 < y < 0.8 and 0.01 < x < 0.1, corresponding to 
a Q 2 range 180 < Q 2 < 7220 GeV 2 . 

Charged particles are selected over a defined range of pt values, where px is 
the transverse momentum with respect to the j*P axis in the HCM system. 
Results for the asymmetry are presented as a function of p c , the minimum value 
of pt of the selected charged particles, since the QCD predictions depend on 
this variable. The measurements are compared to theoretical expectations. 

Previous studies of single hadron production in neutral current DIS have ob- 
served a cos0 term that was attributed to non-perturbative effects J|,||]. Here, 
for the first time, the observation of both cosc/> and cos 20 terms are reported. 
As discussed above, the cos 2<j> term can be unambiguously attributed to per- 
turbative QCD processes. 



2 Experimental Setup and Event Selection 

The data presented here were taken at the positron-proton collider HERA using 
the ZEUS detector. The 38 pb _1 of data were taken in 1996 and 1997, when 
HERA operated with positrons of energy E e = 27.5 GeV and protons with 
energy 820 GeV. A detailed description of the ZEUS detector can be found 
elsewhere |o[0. 

Throughout this paper the standard ZEUS right-handed coordinate system is 
used, in which X = Y = Z = § is the nominal interaction point, the positive 
Z-axis points in the direction of the proton beam (referred to as the forward 
direction) and the X-axis is horizontal, pointing towards the centre of HERA. 
The polar angle, 8, is defined with respect to the positive Z direction. 

The key component for this analysis is the central tracking detector (CTD) |l2j 
which operates in a magnetic field of 1.43 T provided by a thin superconducting 
solenoid. The CTD is a drift chamber consisting of 72 cylindrical layers, covering 
the region 15° < 9 < 164°. The transverse momentum resolution for full- 
length tracks is a(p T )/p T = 0.0058p r © 0.0065 ©0.0014/p T (p T in GeV). The 
interaction vertex is measured using the CTD with a typical resolution along 
(transverse to) the beam direction of 0.4 (0.1) cm. 

Also important for this analysis is the uranium-scintillator sampling calorimeter 
(CAL) [^3|, which surrounds the solenoid. The CAL consists of 5918 cells, or- 
ganised in three sections, the forward (FCAL), barrel (BCAL) and rear (RCAL), 
with longitudinal segmentation into electromagnetic and hadronic sections. The 
energy resolutions, as measured in test beams, are a/E = 0.18/ yj E(GeV) and 
0.35/-\/£ , (GeV) for electrons and hadrons, respectively. 

The angular coverage of the ZEUS detector allows the kinematic variables x and 
y to be reconstructed in a variety of ways using combinations of energies and 
angles of the positron and the hadronic system. Variables calculated only from 
the measurements of the energy, E' e , and angle, 9 e , of the scattered positron 
are denoted with the subscript e, whilst those calculated from the hadronic 
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system measurements, using the Jacquet-Blondel method |14[], are denoted by 
the subscript JB. Variables calculated by these approaches are used only in the 
event selection. In the double angle method p5[ , denoted by DA, the kinematic 
variables are determined using 6 e and the angle jh (which is the direction of 
the struck quark in the QPM), evaluated using the CAL cells corresponding to 
the hadronic final state. 

The DA method was used throughout this analysis for the calculation of the 
boosts and the kinematic variables because it is less sensitive than other methods 
to systematic uncertainties in the CAL energy measurement. 

Corrections are applied |l6| for hadronic energy loss in inactive material in 
front of the calorimeter and for the backscattering of particles into the BCAL 
or RCAL from hadron jets in the FCAL direction. 

The triggering and online event selections were identical to those used for the 
ZEUS measurement of the structure function F% [ p7[ . 

Further selection criteria were applied both to ensure accurate reconstruction of 
the kinematic variables and to increase the purity of the sample by eliminating 
background from photoproduction processes: 

• E' e > 10 GeV, to achieve a high purity sample of DIS events; 

• De < 0.95, to reduce the photoproduction background; 

• Vjb > 0.04, to give sufficient accuracy for DA reconstruction of Q 2 and 

x; 

• 40 < ^2 (E — pz) < 60 GeV, where the summation is over all calorimeter 
cells, to remove photoproduction events and events with large radiative 
corrections; 

• \X\ > 16 cm or \Y\ > 16 cm, where X and Y are the impact position 
of the positron on the RCAL, to avoid the region adjacent to the rear 
beampipe; 

• —40 < Z veltex < 50 cm, to reduce background events from processes other 
than ep collisions; 

• 0.01 < xda < 0.1 and 0.2 < y^A < 0.8, to ensure 7# lies within the CTD 
acceptance. 

The reconstructed tracks used in the charged particle analysis are associated 
with the primary event vertex and are required to have p!f b > 150 MeV and 
|?/ ab | < 1.75, where ?y lab is the pseudorapidity, given by — ln(tan(0/2)). These 
cuts select the region of CTD acceptance where the detector response and sys- 
tematics are best understood. 

Since the gluon fragmentation function is 'softer' than that of the quarks, the 
hadron-quark correlation can be enhanced by selecting 'leading' charged par- 
ticles. This was accomplished by cutting on the Lorentz-invariant variable 
Zh = P ■ Phi P • <7, where p^ is the track or particle four-vector. The tracks 
were selected if 0.2 < Zh < 1.0; no selection biases arise from the r\ cut applied 
in the laboratory frame. 

The number of events that pass the event selection is just over 13,800; approxi- 
mately 7,700 charged tracks satisfy the Zh criteria. In principle, more than one 
entry per event is permitted. However, as can be seen from the above numbers, 
such an occurrence is rare. 
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3 QCD models and event simulation 



Monte Carlo event simulation is used to correct for acceptance and resolution 
effects. The detector simulation is based on the GEANT 3.13 program p8[ . 

To calculate correction factors, neutral current DIS events were generated via 
DJANG06 2.4 @, using HERACLES 4.5.2 §§, which incorporates first-order 
electroweak corrections. The parton cascade was modelled with the colour dipole 
model (CDM), using the ARIADNE 4.10 |2l) program, in which colour coher- 
ence effects are implicitly included in the formalism of the parton cascade. The 
Lund string fragmentation model p2| is used for the hadronisation, as imple- 
mented in JETSET 7.4 ^3|. For the analysis of the data, a Monte Carlo sample 
of 120k events, corresponding to an integrated luminosity of 67 pb _1 , was gen- 
erated with 0.13 < y < 1.0 and 0.007 < x < 0.3 and the CTEQ4D @ parame- 
terisation of the proton's parton distribution functions. The CTEQ4D parame- 
terisation has been shown to describe satisfactorily the HERA measurements of 



the proton structure function Fi 25, 26[ in the kinematic regions investigated. 



For systematic checks, an additional sample of 120k events was generated using 
the DJANG06 and LEPTO 6.5@ Monte Carlo programs. The parton shower 
option, matched to the 0(a s ) matrix elements (ME+PS), was used to simulate 
the parton cascade. The coherence effects in the final-state cascade are included 
by angular ordering of successive parton emissions. For the hadronisation the 
Lund string model was used and for the parameterisation of the parton distri- 
bution functions CTEQ4D was used. The LEPTO sample was generated with 



the soft colour interactions |29| turned off. 



The ARIADNE sample was generated with the high Q 2 modification devel- 
oped during the 1998 HERA Monte Carlo workshop Q. A version of ARI- 
ADNE 4.10, with a corrected implementation of the (j> dependence, was made 
available by the author ]30[ | and it is this version that is used in the comparison 
with the data. The Pomeron-scattering option of ARIADNE was switched off. 



4 Correction Procedure 



The Monte Carlo samples described above were used to correct the differential 
cf> distributions for detector effects. To minimise the sensitivity of the correction 
procedure to the underlying theoretical input of the Monte Carlo generator, the 
data were fitted using the following approach. Rather than correcting the data 
directly, a theoretical distribution T(4>t; a) was used to represent the corrected 
azimuthal angle, <jyp, distribution. The parametrisation T(0^;a) used in the 
analysis is that given in equation |l]) plus an additional Dsmcj) term, so that 
a here denotes the set of parameters A, B, C and D. The set of parameters, o, 
was determined by minimising the quantity 



2 _ (-Pmeas(</ ) M) — Prcc{<t>M\a)Y 



E 



X - / J 2 
U 



M 



where P m ea&{4>M) is the measured cf> distribution and Precis m ] S) is the theo- 
retical distribution smeared for detector effects. The term a 2 M is the sum in 
quadrature of the uncertainties of P mea s and P rec . The detector smearing takes 
account of: (i) migration between different (x, y) bins, (ii) track reconstruction 
efficiency and the 4>t - 4*m migration matrix and (iii) correction factors for 
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pair conversion, K° decay etc. These smearing coefficients were obtained from 
generated and reconstructed Monte Carlo events. 

For a bin size of tt/5 radians, the efficiency for a generated 4> to be reconstructed 
in the same bin as generated is > 55%. The purity of the reconstructed <p bin 
(i.e. the percentage of the tracks reconstructed in a bin coming from particles 
generated in the same bin) is > 60%. 

It should be noted that in this procedure a corrected experimental <j) distribution 
is not determined. Instead the parameterisation of the true T(<f>T] o) distribution 
is fitted directly to the data. An important feature of this method is the fact 
that the (f>T - 4>m migration matrix depends only weakly on the generated (f> 
distribution, and is essentially determined by the experimental inefficiencies, 
which are well simulated. Because significant off-diagonal elements are present 
in the matrix, this fitting technique was used for the main analysis. 

The more conventional bin-by-bin correction method was used as a check. Though 
this method is expected to be more dependent on the underlying physics in the 
generated distribution it has the advantage of producing a corrected experi- 
mental <fi distribution. Correction factors, F(<f>), were obtained from the Monte 
Carlo simulation by comparing the generated distributions before the detector 
and trigger simulation with the "observed" distributions after these simulations 
followed by the same reconstruction, selection and analysis as the real data. 
The correction factor for </> is defined as the ratio of the generated and observed 
differential (f> distributions: 

. . / 1 dn\ I ( 1 dn\ 

where N is the number of Monte Carlo events and n the number of charged 
particles in the (x,y) range. The correction factors are in the range 0.8-1.2. 



5 Systematic Checks 

The major systematic errors can be divided into three types: uncertainties due 
to event reconstruction and selection; to track selection; and to the modelling 
of the hadronic system. No single systematic uncertainty was larger than the 
statistical error in the mean of either cos (f> or cos 2<p. For both mean values, the 
largest effects, which approached the statistical uncertainties, were associated 
with: the inclusion of tracks not associated with the primary vertex; the use 
of different correction procedures for the energy scattered into the calorimeter 
from the interaction of energetic particles in the forward direction; the use 
of an alternative approach to correct for energy losses due to inactive material 
between the calorimeter and the interaction vertex; and the Zh cut being lowered 
by an amount commensurate with its reconstruction resolution. The mean 
value of cos 20 was also reduced by its statistical uncertainty for p c > 1.0 GeV 
when LEPTO was used to correct the data. The upper and lower uncertainties 
extracted from these systematic checks were added separately in quadrature. 

Three additional cross-checks were made on the measurement that were not 
included in the systematic uncertainties. The analysis was repeated using: (i) 
the bin-by-bin correction; (ii) the measured positron energy in place of that 
calculated from the DA method, to check the sensitivity to the boost, and (iii) 
only events without bremsstrahlung from the Monte Carlo sample. In each case 
the results were consistent, within statistical errors, with those from the original 
analysis. 
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6 Results 



Figure |2] shows the differential <fi distributions of charged hadrons for four values 
of the pt cut, p c , in the hadronic centre-of-mass. The results obtained from the 
main analysis method (full line) and the bin-by-bin correction method (points) 
are seen to agree. At low p c , a clear cos0 term is observed. As the value of p c 
is increased a cos 2(f) term becomes evident. 

Figure |^ shows the moments (cos0) and (cos 20) as a function of p c . They 
were calculated from the fit using (cos0) = B/2A and (cos 20) = C/2A. The 
sin0 term is consistent with zero independent of the value chosen for p c . The 
statistical errors are those obtained from fitting the data. The value of (cos <fi) 
is negative and decreases in magnitude as p c is increased. In contrast, the 
value of (cos 2(f>) is positive and rises as p c is increased. Figure || also shows 
predictions of the Lund Monte Carlo programs (LEPTO and ARIADNE). The 
data have larger absolute values of (cos</>) than predicted by LEPTO, though 
for higher values of p c the agreement is reasonable. The dependence of (cos 2(f)) 
on p c is weaker in the LEPTO generator than in the data. Imposing a flat </> 
dependence on the matrix elements implemented in the LEPTO program, in 
order to investigate any asymmetry arising from purely QCD coherence effects, 
results in behaviour different from the data; both moments are almost zero for 
all values of p c . The ARIADNE Monte Carlo predictions tend to have smaller 
absolute values of (cos </>) than the data, whilst the (cos 2(f)) term is in reasonable 
agreement with the data. 

Figure ^ compares the data with two LO QCD calculations. Both calculations 
were made with Q as the appropriate scale, with the Binnewies et al. LO frag- 
mentation function and with the CTEQ4 LO proton parton densities [^4[ . 
The LO calculations result in a qualitatively similar behaviour to the LEPTO 
and ARIADNE Monte Carlo generator predictions. 

The analytic calculation from ZEUS (based on the calculation of Chay et al. ||) 
includes an estimation of the non-perturbative contribution, from intrinsic kx 
and hadronisation pt, and integrates over the whole kinematic range. The 
results of Ahmed & Gehrmann are purely perturbative at leading order in a s 
and are evaluated at the mean values (x) — 0.022 and (Q 2 ) — 750 GeV 2 of the 
data. The different implementations account for the observed difference in the 
two predictions; using (x) and (Q 2 ) in the ZEUS perturbative calculation leads 
to agreement with the Ahmed & Gehrmann calculation. 

There are perturbative and non-perturbative contributions to both (cos (f>) and 
(cos 2(f)). For (cos0), the non-perturbative contribution to the ZEUS calcula- 
tion is at most 20% at low p c for a mean intrinsic kx of 0.5 GeV and a non- 
perturbative fragmentation px of 0.5 GeV. However, this fraction depends on 
the values of kx and px used. For (cos 2(f)) , this non-perturbative contribution 
is negligible, even for mean values of kx = 0.9 GeV and px — 0.9 GeV. It can 
be concluded that the observation of a significant (cos 2(f)) term is clear evidence 
for a perturbative contribution to the azimuthal asymmetry. 

The sensitivity to higher-order corrections was examined at the partonic level 
using the DISENT Monte Carlo program (3|] . This study suggests that higher- 
order corrections are large (factors of 2-3) at lowp c for the (cos (f>) term, resulting 
in larger negative values. For the (cos 2(f)) term they are much smaller. This 
may explain the disagreement between the data and the LO QCD models for 
(cos (f>) at low pc. 
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7 Conclusions 



The azimuthal asymmetries in the deep inelastic electro-production of single 
particles have been measured at HERA in the hadronic centre-of-mass frame in 
the kinematic region, (x) = 0.022 and (Q 2 ) = 750 GcV 2 . For hadrons produced 
at large transverse momenta, the measured value for (cos <fi) is negative and is 
in agreement with QCD predictions. The moment (cos 20) has been measured 
here for the first time and is non-zero and positive. It increases as a function 
of the minimum particle transverse momentum, as expected from QCD. Since 
the non-perturbativc contribution to (cos 2<p) is predicted to be negligible, this 
measurement provides clear evidence for a perturbative QCD contribution to 
the azimuthal asymmetry. 
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Figure 1: Definition of the azimuthal angle (p. The incoming lepton is denoted 
by e, the scattered lepton by e', the exchanged boson by 7* and the outgoing 
hadrons or partons by h. The dotted line represents the intersection of the e — e' 
scattering plane with the transverse plane. 
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Figure 2: The differential distributions obtained for four values of the pr cut, 
p c , in the hadronic centre-of-mass frame in the kinematic region 0.01 < x < 0.1 
and 0.2 < y < 0.8 for charged hadrons with 0.2 < zh < 1.0. The full line 
(with accompanying statistical error band) is the result obtained from the main 
analysis using the unfolding technique. The data points were corrected using a 
bin-by-bin procedure (only statistical errors are shown). 
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Figure 3: The values of (cos</>) and (cos2</>) are shown as a function of p c in 
the kinematic region 0.01 < x < 0.1 and 0.2 < y < 0.8 for charged hadrons 
with 0.2 < Zh < 1.0. The inner error bars represent the statistical errors, the 
outer are statistical and systematic errors added in quadrature. The lines are 
the predictions from LEPTO (full line) , LEPTO with flat 4> dependence (dashed 
line) and ARIADNE (dotted line) Monte Carlo models. 
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Figure 4: The values of (cos </>) and (cos 2</>) are shown as a function of p c in 
the kinematic region 0.01 < x < 0.1 and 0.2 < y < 0.8 for charged hadrons with 
0.2 < Zh < 1.0. The inner error bars represent the statistical errors, the outer 
are statistical and systematic errors added in quadrature. The lines are the LO 
predictions from ZEUS with perturbative and non-perturbative contributions 
(full line), ZEUS with the perturbative contribution only (dashed line) and 
Ahmed & Gehrmann (dotted line - see text for discussion). For the case of 
(cos 20), the ZEUS total and perturbative predictions are almost identical. 
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